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Abstract 

A sequence of constant mean curvature surfaces Ej with mean curvature Hj — > cxd in a three- 
dimensional manifold M condenses to a compact and connected graph T consisting of a finite 
union of curves if Sj is contained in a tubular neighbourhood of T of size 0{1/Hj) for every 
j S N. This paper gives sufhcient conditions on T for the existence of a sequence of compact, 
embedded constant mean curvature surfaces condensing to T. The conditions are: each curve 
in 7 is a critical point of a functional involving the scalar curvature of M along 7; and each 
curve must satisfy certain regularity, non-degeneracy and boundary conditions. When these 
conditions are satisfied, the surfaces Sj can be constructed by gluing together small spheres of 
radius 2/Hj positioned end-to-end along the edges of T. 

1 Introduction 

Background. A constant mean curvature (CMC) hypersurface S contained in an ambient (for 
simplicity three-dimensional) compact Riemannian manifold M has the property that its mean 
curvature with respect to the induced metric is constant. This property ensures that area of S is a 
critical value of the area functional for surfaces of M subject to an enclosed- volume constraint. The 
study of CMC surfaces in the classical setting of M = is a well established field of Riemannian 
geometry and the literature concerning the construction and properties of such surfaces is enormous. 
Indeed, it is known that there is a very great flexibility in the construction of CMC surfaces and 
that also a number of interesting rigidity results hold as well. 

Examples of such rigidity results are as follows. In the particular case of finite-topology sur- 
faces, Meeks [11] proved that any end of a complete Alexandrov-embedded CMC surface in 
is cylindrically bounded; while Korevaar, Kusner and Solomon [8] improved this by showing that 
any such end converges exponentially to one end of a Delaunay surface. Furthermore, the possible 
directions of the axes of these ends are also subject to limitations, as well as the flexibility to change 
these directions within the moduli space of all such surfaces, see [21 [3l H]. These limitations are 
phrased in terms of a certain flux that was discovered by Kusner. The flux is a vector associated 
to any closed loop in a CMC surface; it is constant under deformations of this loop, and in fact 
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only depends on the homology class of this loop in the surface. There is a flux associated to a 
simple positively oriented loop around each asymptotically Delaunay end, which depends only the 
direction of the axis and the neck-size of the limiting Delaunay surface. The homological invari- 
ance also shows that the sum over all ends of these limiting fluxes must vanish, which is a global 
balancing condition for the entire CMC surface. One immediate consequence is the non-existence 
of a complete Alexandrov-embedded CMC surface in with only one end. Finally, Korevaar and 
Kusner [7] have a very general structure theorem for these types of CMC surfaces which states 
that the ends are as described above while the remainder of the surface is contained in a tubular 
neighbourhood of a graph consisting of line segments and a balancing formula holds at the nodes 
of the graph. 

The corresponding picture in general Riemannian manifolds is considerably less well-developed. 
Rather than focus on special cases (such as when M is the sphere or hyperbolic space, where there 
are many interesting results) we will work in an unspecified Riemannian manifold but consider small 
CMC surfaces of high mean curvature. In this setting, Rosenberg |13j has shown that if S is a 
closed CMC surface in an arbitrary (compact) 3-manifold M, with sufficiently large mean curvature 
H, then M \ S has two components, and the inradius at any point in one of these components is 
bounded above by C/H. In other words, S is localized in a small tubular neighbourhood of T. 
This phenomenon is captured via the following definition. 

Definition 1. A sequence of constant mean curvature surfaces in M with mean curvature 
Hj —5- cxD condenses to F if Ej is contained in a tubular neighbourhood of F of size 0{1/Hj) for 
every j G N. 

We are led to the following two central and very natural questions. 

• What are the possible condensation sets F for condensing sequences of CMC surfaces Sj? 

• If condenses to F then what can be said about the nature of for large jl 

As an example of the kind of answers we hope to provide, let us consider a sequence of dilations of 
a complete CMC surface with k asymptotically Delaunay ends in M'^, where the dilation parameter 
goes to zero. Such a sequence condenses onto a one-dimensional set, here a union of half-lines 
meeting at a point where the unit vectors defining the half-lines satisfy the fiux balancing condition 
described above (which specifically means that the sum of these unit vectors, weighted by a factor 
relating to the Delaunay parameter of the corresponding end, vanishes). 

An obvious hope is that this Euclidean example represents the local behaviour of condensing 
sequences of CMC surfaces in general Riemannian manifold; i.e. that in a Riemannian manifold M, 
a condensation set F is some sort of network of geodesies, where each edge of this geodesic network 
to have a 'weighting' which carries information about the Delaunay parameters of the CMC tubular 
piece which converges to that edge, and that the nodes of this network are balanced. There is in 
fact some evidence that this is true, in the form of a result of Mazzeo and Pacard [9]. They show 



2 



that if r is any closed, non-degenerate geodesic in M, then most geodesic tubes of sufficiently small 
radius about T can be perturbed into CMC surfaces with large H (thus providing an example of 
a sequence of CMC surfaces condensing to T. Furthermore, they show that if Sj condenses onto a 
smooth one-dimensional manifold T and the pointwise ratio of the norm of the second fundamental 
form to the mean curvature is bounded above independently of j, then T must be a geodesic. 

Despite this compelling evidence, however, it seems that the actual characterization of condens- 
ing sequences of CMC surfaces in M is potentially far more complex and interesting. To understand 
the reason why this is so, it is helpful to review the first result about CMC surfaces of high mean 
curvature. This result is due to Ye [H] in the early 1990s and relates to sequences of CMC surfaces 
collapsing to a single point, so that F = {p} for some p in M. Ye's result has two components. First, 
he proves that if p is a non-degenerate critical point of the scalar curvature of M, then geodesic 
spheres around p with small radius may be perturbed into CMC surfaces with large H. Second, 
he also proves conversely that if a sequence of CMC spheres with H ^ oo forms a local foliation 
with uniformly bounded pointwise ratio of the norm of the second fundamental form to the mean 
curvature, then this sequence converges to a critical point of the scalar curvature. 

The conclusion to be drawn from Ye's work is that the scalar curvature of M may be involved 
in the behaviour of condensing sequences of CMC surfaces in some way. The question now becomes 
if the scalar curvature of M influences the behaviour of sequences of CMC surfaces condensing to a 
one-dimensional variety. Given Mazzeo and Pacard's result [9] cited above, we should only expect 
such an influence if F is a non-smooth variety or else the pointwise ratio of the norm of the second 
fundamental form to the mean curvature of Sj diverges on at least one sequence of points pj € Sj. 

Butscher and Mazzeo [T] have recently shown that it is indeed possible to construct a condensing 
sequence of CMC surfaces where the condensation set is contrary to what one would expect in the 
Euclidean case. Namely, they construct surfaces that locally resemble a Delaunay surface but are 
either compact or have one end. It should be noted that their construction depends critically on 
the nature of M — they assume that M is cylindrically symmetric (about a geodesic 7) and very 
rapidly asymptotically Euclidean. They use the gluing technique to construct their surfaces: a 
collection of 0{L/r) small geodesic spheres of radius r are positioned end-to-end along 7 and glued 
to each other using re-scaled, embedded catenoids to produce a compact, almost-CMC surface of 
length 0{L). A semi- infinite embedded Delaunay surface can also be glued to one of the terminal 
spheres of this configuration to produce a one-ended almost-CMC surface. These surfaces are 
then perturbed into true CMC surfaces by solving a partial differential equation. The sizes of the 
catenoidal necks employed in the construction vary along the length of 7 in a manner governed by 
a flux formula that implies that the difference between successive neck-sizes can be expressed in 
terms of the gradient of the scalar curvature along the axis connecting these two necks. Moreover, 

^The gluing technique is a well-known technique used in geometric analysis. It goes back to the work of Kapouleas 
[51 |B] in the context of CMC surfaces and has been further developed by many other researchers. See |10l 112] for 
surveys about the current state of this approach 
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these neck sizes are quite smah compared to the radii of the spheres — the pointwise ratio of the 
norm of the second fundamental form to the mean curvature diverges precisely in the neck regions 
of these surfaces. 

The condensation sets in Butscher and Mazzeo's work above are very simple: a geodesic segment 
on the one hand and a geodesic ray on the other. The reason for this is the highly restrictive 
symmetry conditions satisfied by M, which in particular forces the gradient of the scalar curvature 
to point along 7. One might thus hope for more interesting behaviour and more radical departures 
from the Euclidean picture if these symmetry conditions are weakened. In this paper, we consider 
a completely general (though compact) ambient manifold M containing a union of curves F. We 
identify a set of sufficient conditions on the curves in T and on the endpoints of these curves that 
must hold in order for the construction of a sequence of embedded CMC surfaces condensing to T 
to be possible via the gluing technique. These conditions, expressed in greater detail below, show 
that the scalar curvature affects the curves in T as well as the balancing formulae satisfied at the 
endpoints of these curves. In particular, these curves are in general neither geodesies nor integral 
curves of the gradient of the scalar curvature. 

Results. Consider a graph F in M with edges £ := {71, . . . ,7s} and vertices V := {pi, . . . ,Pn}- 
We will assume that the vertices are found amongst the endpoints of the edges and never amongst 
the interior points of the edges. Furthermore, if there are two edges 7e and 7e' emanating from pi 
then the one-sided tangent vectors of 7e and 7e' at pi are not co- linear. We will also assume that 
each 7e is parametrized by arc- length and we will also use the mapping t 1— > 7e(i) to denote this 
parametrization. We now introduce four existence conditions on the curves in F. 

1. Each 7e G <S is smooth with geodesic curvature uniformly bounded independently of r. 

2. For each 7 € f and all sufficiently small r, there is a smooth, positive function / : 7 ^ M so 
that the pair (7, /) satisfies the equation 



where is a geometric constant and R is the scalar curvature function of M. 

We note here that (jlap is actually equivalent to a system of two equations that we can obtain 
by expanding and projecting parallel and perpendicular to 7. That is. 



f{t) = Qr^{R o 7(t) + c). Hence it is necessary only to solve the first of these two equations. 
3. The following boundary conditions hold for the various pairs (7, /). 





Since (Vflo7,7) = ^ 



/V^7 = nr^{VRojy 

df , (lb) 
dt 

(^R o 7) the second equation implies that there is a constant c so that 
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• For each p E V such that there is only one edge 7 € £" emanating from p, we let / be the 
function corresponding to 7 and assume without loss of generality that 7 is parametrized 
to begin at p. Then 

/(0) = aswehas 7(0) || Vi?(p) and (7(0), Vi?(p)) < . (2a) 

• For each remaining p £ V, we let 71, ... , 7/^ £ he the K > 2 edges emanating from p 
and we assume without loss of generality that these curves are parametrized to begin at 
p. U fi, . . . , Jk are the corresponding functions, then 

K 

(/i(0) + ^7r3(Vi?(p),7^(0)))7^(0) = nr^VRip) . (2b) 

i=l 

4. Let be a continuous one-parameter family of graphs such that — > F where we mean the 
standard notion of convergence for one-dimensional varieties. Then each curve 7t € F^ with 
7( — > 7 € F carries an infinitesimal deformation vector field X^. We will also assume that no 
belongs to the kernel of the linearization of the operator given in (jlap corresponding to 7. 

The main theorem that will be proven in this paper is that a graph F satisfying the conditions 
(1) - (4) above is the condensation set of a sequence of CMC surfaces. 

Theorem. Let T be a graph satisfying conditions (1) - (4) above. Then there is a sequence rj 
and a family of CMC surfaces Ej where each T,j has mean curvature 2/rj and is contained in a 
tubular neighbourhood of radius 2rj ofT. 

The proof of this result can be outlined as follows. We use the standard machinery of gluing 
constructions, which constructs the CMC surfaces condensing onto F in several steps. First, we 
position spherical building blocks (well-chosen small perturbations of spheres) of radius r end-to- 
end along each of the edges in F and glue them to each other by means of optimally matched, 
re-scaled and truncated catenoids. This yields an approximately CMC surface Sr(F). We also 
construct neighbouring approximately CMC surfaces whose spherical constituents and necks are 
displaced by arbitrary small amounts. Next, for each of these approximately CMC surfaces, we 
solve the constant mean curvature equation up to a remainder term in the approximate co-kernel 
of the linearized mean curvature operator of Er(F). We then choose F so that the highest-order 
term (in powers of r) of the remainder vanishes — which is known as finding a balanced initial 
approximately CMC surface. Finally, we vary this F slightly in order to find a surface for which 
the remainder term vanishes identically. 

The key feature of the proof above is that Conditions (1) - (4) given above translate into the 
equations that must be solved in order to find the graph F that yields a balanced initial approx- 
imately CMC surface. To explain this in more detail, we must first describe how the balancing 
equations are derived. First let us consider a CMC surface S with mean curvature h in (M, g) and 
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suppose that U and W are open sets in S and M, respectively, such that d\V = Z7 U Q for some 
surface- with-boundary Q. The first variation formula for the area olU with the volume of W fixed 
relative to the one-parameter family of diffeomorphisms (j)t generated by a vector field V on M then 
implies 



0= ^^[Area[4>t{^nU)) -hVol{(t)t{W)) 



g{u,V) + h / g{N,V), (3) 

t=0 JdU 



where v is the unit normal vector field of V in S and N is the unit normal vector field of Q in M. 
Next we consider the case when S has approximately constant mean curvature, as is the case for 
our approximate solutions Sr(r). It turns out that the perturbation argument for finding a nearby 
exactly CMC surface works only when the right hand side of ([3]) is sufficiently close to zero for 
all vector field V in the approximate co-kernel of the linearized mean curvature operator. Since 
Sr(r) for sufficiently small r consists of spherical pieces joined by small necks, these vector fields 
are the approximate translation vector fields (defined as the coordinate vector fields of a normal 
coordinate chart centered at a spherical piece) multiplied by cut off functions that cause them to 
vanish on all but one of the spherical pieces. Evaluating the right hand side of ([3]) on such a vector 
field corresponding to the spherical piece centered at p G F yields 

r.h.s.= r£i{r]i,V) -Qr^{VR{p),V) + 0{re^) + 0{r^) (4) 

all necks 

where r/j is the unit vector pointing along the geodesic from the center of the spherical piece in 
question to the i^^ neck that connects this spherical piece to its neighbours, rsi is the width of this 
neck and e := maxj{ej}. Also, ^R{p) is the gradient of the scalar curvature of M at p, the brackets 
(•, •) denote the ambient metric at p and is a constant independent of r. 

In order to find F so that the right hand side of ^ vanishes identically, we proceed as follows. 
Let V denote the set of centers of all the spherical pieces used in the construction of Sr(r)- We 
first find F so that the leading order terms in Q all vanish, namely that the equations 

0= ^ e^r^i-nr'^VRip) (5) 

all necks 

hold for all p € V. We then perturb F slightly so that the remaining small errors can be made 
to vanish as well. This latter step requires a non- degeneracy condition: that the map which takes 
graphs near F to values of the right hand side of equation ^ be locally surjective. 

The requirements for the procedure above can be translated directly into Conditions (1) - (4) 
of the Main Theorem. To do so, consider the family of equations of the form ([5]) corresponding to 
the interior of a given curve 7 C F. Dividing through by r, this equation reads 

= i + e'T]-) - ^r'^VR{p) (6) 

where Vj^ are the unit vectors of pointing in the direction of the neck ahead (+) and behind 
(— ) the point p on the curve 7 while are the corresponding neck sizes. We can now view this 
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equation as the discretization of a differential equation. Up to additional small error terms, we have 
r/^ = ±7(p ± r) and if we introduce a neck-size function / : 7 — s- M then we have = f{p it r). 
Hence the right hand side of Q is equal to 

V.,{fj)-nr^VR{p) 

up to additional small error terms. The vanishing to highest order in r of ([6]) is thus equivalent to 
Condition (2) of the Main Theorem. Condition (1) is now the regularity condition which ensures 
that the various estimates of 5]r(r) that are needed in the perturbation theory remain uniform. 
Condition (3) is obtained by expressing the equations that must hold at the boundaries of all the 
curves in T as boundary conditions for the the pairs (7,/)- Finally, Condition (4) is the required 
non-degeneracy condition. 

Comments on the existence conditions. The task of finding a non-trivial network of curves 
r satisfying the existence conditions above is of course still open. (The curve used in Butscher- 
Mazzeo [1] should be seen as a trivial solution of these conditions because there the gradient of 
the scalar curvature points along the curve.) We will not attempt to solve in any general way the 
existence conditions in this paper. Instead, we will point out a key feature of the equation (jlap 
which should serve as a starting point for the investigation of the existence and properties of its 
solutions. 

Proposition 2. Curves satisfying equation (jlap are critical points of the functional 

7f + 17V(iio7 + c)2) 
for some constant c. 

Proof. We multiply both sides of (fTa|) by / and re-parametrize 7 so that fj:=& and hence 
\\a\\ = f = Qr'^{R o 7 + c). This yields the equation 

V^d = nr^\\&\\VRo^ = OV^(i?o7 + c)Vii 07 = 17V^V(ii o 7 + c)^ 

which is precisely the Euler-Lagrange equation of the given functional. □ 

Acknowledgements. I would like to thank Rafe Mazzeo, Frank Pacard, Harold Rosenberg and 
Richard Schoen for interesting discussions during the course of this work. I would also like to thank 
Harold Rosenberg for his generous hospitality during my stay at IMPA in Brazil in May 2009. 

2 Constructing a Family of Initial Surfaces 

Let r be a graph as defined in Section [1] with edges £ := {71, . . . , 7^;} and vertices V := {pi, . . . ,Pn}- 
We now show how to construct a family of initial surfaces based on F by positioning small slightly 
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deformed geodesic spheres of radius r end-to-end along the edges in £, and gluing these to each 
other by inserting small catenoids. Once a member of this family of surfaces is chosen, we also show 
how to construct a family of perturbations of this surface by allowing the locations of its spherical 
constituents to vary. The construction procedure can be explained in several steps. 

Step 1: The building blocks. We first construct building block surfaces in W^. Identify 
with the unit sphere centered at the origin in and let Jg := x^\^2 be the restriction of the s*'* 
coordinate function to Set }Cg2 := span{Ji, J2, J3}. Choose a collection of points qi, . . . ,qn 
and small parameters ei,...,e„ G M+ (which we'll abbreviate here by q and e) and define the 
function G^^g : \ {qi, . . . , qn} — > M as the unique solution of the distributional equation 

n 

As2 G^^g + 2G^^g =^^ei5i + J and G^^g _L ^2 /C§2 (7) 

i=l 

where 5i is the Dirac 5-mass at pi , while Ag2 is the Laplacian of with respect to the induced metric 
and J G /C is such that the right hand side of d?]) is L^-orthogonal to /Cg2, thereby guaranteeing 
the existence of G. Now let e = maxjei, . . . and let := e^/^ (the reason for this choice will 
become clear in Step 2 below). Define 

Sr[qA ■■= + G,>-(0))0 : G §2 \ IJ i?,^(g,)| 

which is the normal graph over \ IJiLi ^r^{qi) generated by the function G^^g. 

We now transplant these building blocks to M. Choose a point p ^ M and an orthonormal 
frame E := i?2; -E's} at P from which we construct the inverse of the geodesic normal coordinate 
map (t)~\ : ^ M via by (t)~\{x) := expp(x^£^i x^E^) for x = (x^, x^, x^) G -B_r(0) where R 
is some radius smaller than the injectivity radius of M at p. This allows us to obtain a building 
block in M via the prescription 

Sr^,E,q,e\ := (t)~^j^{Sr[q, e]) . 
Finally, let us refer to the function G^^g as the generating function of the building block Sr\p,E,q,£]. 

Step 2: Gluing two building blocks together. Let us suppose that we are given one of the 
building blocks defined above, say S := Sr[p, E,(j,e\ with q = (gi,...,g„) and e = (ei,...,e„). 
We now show how this building block can be glued to a second building block at the boundary 
component of S corresponding to a chosen q^. First, view qi as a unit vector in TpM by means of 
the identification q^ := Y2'j=i Qi^j and parametrize the geodesic emanating at p in the direction Qi 
by o"i(t) := expp{tqi). Choose the point p' := ai{{2 + T)r) where r > is a parameter satisfying 
r = 0{elog{l/e)). Next, choose an orthonormal frame E' := {E[, E'2, E'^} at p' such that E[ = 
-&i{{2 + T)r). Finally, let S' := Sr[p' , E' , q' , e'] where q" := {q[, . . . ,q'^,) has q[ := (1,0,0) while 
^2, • • • , g^/ and e' := {e'l, . . . , e'^,) are at this point arbitrary. 
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Letp := (7j((l-|-T/2)r) denote the point half-way between p and p' and let F := &i{{l+T/2)r). 
Choose an orthonormal frame at of the form := {F'', .Eg, £^3} and consider the image of a 
neighbourhood of that contains both S and S' under the geodesic normal coordinate mapping 
(f)p\, gb. Re-scaling the image of this neighbourhood by a factor of l/r, we find that p^ is mapped 
to the origin, the geodesic segment (jj is mapped to the x^-axis and the images of S and S', near 
the origin when r is small, can be represented as two graphs over the (x^, a;^)-plane of the form 

{{x\x^x^):x^ = F,,i,{\\{x^x'm and x^, x^) : = F^JIKx^, x^)]])} 

respectively. One can check that the Taylor series expansions of the generating functions of S and 
iS" near their singular points imply expansions for Fsp^ and F^^^ near zero of the form 

Kphix) = 1 + ^1 (c' + C log(x)) + 0(x2) + 0{\e'\x'') 
Fsph{x) = - ei(c + Clog(x)) + 0{x^) + 0{ex^) 

for X > 0. Here, c, c/, C, C are constants. 

We will now show how to construct a six-parameter family of interpolating surfaces between 
the images of S and S' by gluing S and S' together using rotations, translations and re-scalings of 
a standard catenoid (i.e. the cylindrically symmetric, two-ended minimal surface in R^). We first 
show how to construct the "optimal" member of this family. To begin, consider a catenoid with its 
axis of symmetry along the x^-axis, scaled by a factor of , and translated by an amount d' along 
this axis, where and d}' are yet to be determined. The x^ > end of this catenoid is given by 

= ^nef'^(^''^') •= e^arccosh(||(x^x=^)||/e^) + 

= e^(log(2)-log(e^))+e^log(||(x^x=^)||)+£^c^^ + 0((e^)Vll(x^x3)||2) 

near the origin; and the x^ < end of this catenoid is given by 

= ■■= -£^arccosh(||(x2,x3)||/e^) +e^d^ 

= -£^(log(2)-log(e^)) -£^log(||(x2,x3)||)+e^d^ + 0((£^)Vl|xf) 

near the origin. Optimal matching of the asymptotic expansions of the ends of the catenoid with 
the asymptotic expansions of the images of S and S' then requires that the equations 

hold, where 

k{e') := (2(log(2) - log(e^)) " ^ " ^) • (9) 

These equations imply that and r together determine the parameters of the optimally matched 
neck and as well e'^. In other words, the building block S and the spacing r determines the 
parameters of the neighbouring building block S' at the point of connection with S as well as the 
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neck interpolating between S and 5". Note also that the matching between the images of S and S' 
and the neck defined by this choice of parameters is most optimal in the region of the (x^, x^)-plane 
where the error quantity + 0((e^)^/||(x^, is smallest. It is easy to check that 

this occurs when = 0(max(e'^/^, (e')^/^)). To complete the gluing process, we define a 

smooth, monotone cut-off function x ■ [0, oo) — > [0, 1] which equals one in [0, ^] and vanishes outside 
[0, 1]. Next, we define the functions : Bi{0) ^ M by 

F^{x^x^):=x{\\{x^x^)\\/{e'f/')Fif^'^{x^ (10) 

where F*^^ = Fsph if * is — and F*^f^ = F^'^^ if * is +. Now define the interpolating surface as 

N:=cP^,]^,[{{F-^{x^x'),x\x^):e' < \\{x\x^)\\ < {e'f/^} 

U{{F-{x\x^),x'^,x^) : < \\{x\x^)\\ < (e^)^/^}) . 

The surface obtained by gluing S to S' can therefore now be defined as S U U S". 

Next, we explain how to construct interpolating surfaces close to the optimal one we have just 
defined, where the choice of catenoid neck can vary in a six-parameter family of choices. Let 
A'^ denote the unit-scale catenoid with its axis of symmetry along the x^ axis and centered on 
the origin. Then the catenoid used in the optimal interpolation is thus (^N + ((i'',0,0)). Let 
Re2fiz be the rotation by an angle 62 about the x^-axis followed by a rotation by angle ^3 about 
the x^-axis. Let e € M be a small real number and let d := (^1,^2,^3) € be a vector of 
small norm. Define the catenoid N[di,d2,dz,02^6^,e\ := (1 -|- e)e^ Re2fiz{N + (d'',0, 0) + d^ If 
\s\ + 1^21 + l^sl + \\d\\ is sufficiently small, then the x-^ > and x^ < ends of N[di,d2-,d^,e,62.,9^] 
can still be written as graphs over the (x^, x^)-plane. We then define A^[(ii, (i2, (^3, 6*2, 6*3] as the 
interpolating surface obtained by replacing the functions -F^g^j. ' by the graphing functions of 
N[di,d2-,d-i,62.,9^,e\ in the definition of equation ([TO]) . We get our six-parameter family of glued 
surfaces S U N[di,d2-, d^,e, 6*2, 6*3] U S' in this way. 

Step 3: Constructing the initial surface. We first produce a family of piecewise-geodesic 
approximations of T that are parametrized by a collection of small, positive separation parameters 
ri*^ associated to each 7e G 6. We can do this inductively as follows. Let 7e € and ri"^\ . . . , Te^"^ 
be given and let q^"^ = 7e(0) be the vertex at the beginning of 7e. Now determine gi^"* as the unique 
point in 7e located a distance (2 -|- Te^'^)r from q"e^ . Then determine qif^ as the unique point in 7e 
located a distance (2 -|- from (7e^\ Proceed in this way until we reach the vertex at the end of 

7e, namely qe " = 7e(|7e|) where |7e| is the length of 7e. Note that the last separation parameter 
Te^"^ is a function of |7e| and r and t^^'' for s = 1, . . . , 5e — 1 and furthermore, if ri^"^ is to be of 
size o(r) then we must choose r belonging to a sequence of small intervals accumulating at zero 
and a corresponding diverging sequence of integers Se = Se{r). 

We can complete the construction of the initial surface based on the piecewise-geodesic ap- 
proximation found above. We need only position the appropriate building block at each and 
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connect them with the appropriate necks in order to end up with a smooth surface that we will 
call Sf.(r) with r indicating the dependence on the separation parameters chosen above. Indeed, 
if qf^ is attached to neighbours qi by geodesic segments of length (2 + Tj)r for i = 1, . . . , n, then 
we use the building block Sr{q^\qi, ■ ■ ■ ,qn, /Ci, . . . , Kn{Tn) /Cn) where Aj and are the 

quantities appearing in ^ corresponding to qi. The separation parameters also determine the 
optimal interpolating catenoids as we have seen. 

Step 4: Constructing perturbations of the initial surface. Let V := {qi, . . . ^q^} be the 
union of all points in all edges of T where we have placed building blocks during the construction 
of T,r{T) in Steps (1) - (3) above. Note that V C V. For each g G V introduce a tangent vector 
Wq € TqM of small length and denote hy W ^ \[^^^TqM the quantity W := {Wq^,. ..). Form a 
new piecewise-geodesic approximation of F as follows. First move each q to expq(rWq). Then 
for each pair q and g' £ V that are connected by a geodesic segment in the piecewise-geodesic 
approximation of F constructed in Step 3, we connect the corresponding points expg(rPFg) and 
expg/(r/g/) by a geodesic segment. Now we can proceed as in Step (3) and place building blocks at 
each point exp^(rVFg) for Q' € V, where the various parameters of this building block are determined 
by the new direction vectors and the lengths of the geodesic segments attached to expq(rWg). For 
each G V , introduce H^b G of small length and denote by H € n^^^— b the quantity 
S := (S^b,...). Now glue the building blocks adjoining the edge containing together using 
the interpolating catenoid whose parameters are given by the components of H^b. Once we've 
completed the process that we have just described, we end up with a smooth surface that is a small 
perturbation of the surface constructed in Step 3. Denote this surface by Sr(F, r, W, H). 

3 Deforming an Initial Surface into an Almost-CMC Surface 

Let S.r := Sr(F, r, W, H) be a fixed, initial surface constructed as above. In this section of the paper, 
we show how to deform into an almost- CMC surface by solving a partial differential equation 
using a contraction mapping argument in a space of C^'° functions on S,,- This means that we 
will be able to find a small normal deformation of that is determined by the solution of the 
PDF which yields a surface whose mean curvature is the constant 2/r plus a small, well-controlled 
error term. The arguments required here are all fairly standard; in particular, they are carried out 
in full detail in [1] in the special case where F consists of one curve and the ambient manifold M 
possesses a high degree of symmetry. Hence the arguments will be stated here in an abbreviated 
sense for the convenience of the reader. 

3.1 Function Spaces and Norms 

Partitions of unity. We define once and for all a family of partitions of unity for that will 
be used throughout the paper. In what follows, let V be the set of end-points of all the geodesies 
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segments used in the construction of as above, and also let V be the set of mid-points of all 
these geodesic segments — which is where the various necks have been placed. Let be the scale 
parameter used in the definition of the neck placed at q''' G v'' and put e = maxje^b : q' € V**}. 
We begin by defining subsets of S^- First, for any p e M and < a < ctq for some fixed, small 
threshold (To, we define AnUaip) ■= B^{p) \ B„i2{jp)- Now let g'' G V and define the subsets: 



and 



and write := ^^'^ ^ ^^'^ where ^^'^ the two disjoint components of 'ahead' and 
'behind' the point with respect to the parametrization of the edge to which q^ belongs. If 
o"o is sufficiently small, then the set \ U^bgy^ ^A/J^ U ^ is the disjoint union of open sets 
corresponding to each g G V. In this way we can unambiguously define Sg via the prescription 



Wc now define the smooth, positive cut-off functions making up the partition of unity subordi- 
nate to the cover of given by the sets J^^, and for one fixed o" < ctq. Given (f 
and (/ G V connected to neck regions centered at for i = 1, . . . , n, we define 



Interpolation 




X 



x G 

elsewhere 



Interpolation 




with the property that E^evX^,? + S.^gv^ ^leck,q^ 



X G 

xGU?=irr 

elsewhere 
1 for all a. 



Weighted norms. We introduce the standard weighted norm that achieves control of functions 
in the neck regions and asymptotic ends of S,,. We first define a weight function (^j. : S,, ^ M as 
follows. For each € choose (T{ct) := ©(re^b log(l/|egb |)) on the order of the width of the 



neck regions at . Let 5. 

.a{q') 



neck,q^ 



be such that S^eck,q^{^) gives the distance of a point 



G AA^''* ' to the narrowest part of Ml^'^ We now define the weight function by 



Up) :-- 



4 + i^neck,q^ix)? 



X 



Interpolation 
r 



X G n^''^) 



elsewhere 
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where the interpolation is such that Cr is smooth and monotone in the region of interpolation, and 
has appropriately bounded derivatives. 

Next we define the norms we will use. First, if U is any open subset of and A is any tensor 
field on lA, define 

1^10,,,:= sup II and [A^u := sup "^^""l" p'^'/.f , 

x&u x,x'&u dist(3;,x')" 

where the norms and the distance function that appear are taken with respect to the induced metric 
of Tir, while '^x,x' is the parallel transport operator from x to x' with respect to this metric. Here 
we use the convention that | • |o,w = [']a,u = if ^ = 0. For any function f -.U we now define 
its C^'" norm by 

k 
1=0 

Function spaces. The function spaces that will be used in the remainder of the paper are 
simply the usual spaces C^''^{Er), but endowed with the Cu'"" norm defined above. This space will 
be denoted The norm will often be abbreviated by | • |„fc,Q when there is no cause for 

confusion. 

3.2 Setting Up a Contraction Mapping Problem 

Let ^ : C^'"(Sr) Emb{T,r, M) be the exponential map of in the direction of the unit normal 
vector field of T,r with respect to the backgroung metric g. Hence Hrfij^r) is the scaled normal 
deformation of Sj. generated by / € C^'"(Sr.). The equation 

H[M{^r)]=l (11) 

selects / G C^'"(I]r.) so that firfi^r) has constant mean curvature equal to ^. At this stage, it will 
not be possible to solve (jlip exactly, but using a contraction mapping argument together with the 
choice ly G (1,2) the equation (|lip can be solved up to a co-kernel error term. This means that a 
solution of 

can be found, where £ belongs to a subspace of functions associated to the approximate co-kernel 
of the linearized mean curvature operator of Ti^ that will be denoted /C and specified below. At 
first glance, the error £ will come from terms in the solution procedure that are not sufficiently 
small. However, the true reason for the presence of £ is geometric and will be explained in Section 
m where we show how to eliminate it by a suitable choice of W. 

We now outline how this will be done. First we decompose the mean curvature operator of the 
deformed surface as 

H[firf{tr)] = H[tr]+C{rf) + Q{rf)+n{rf) 
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where £ := A + \\B\\'^ is the dominant part of the hnearized mean curvatm'e operator, Q is the 
quadratic remainder part of the mean curvature and 7Y(r/) := Ric\^ ^g, )(^/' ^f) additional 
smah error term. Then we construct a bounded parametrix TZ : C'^'^2{^r) Cu'^C^r) satisfying 
C o TZ = Id + £ where £ maps into )C. Now the Ansatz 

f:=lnl^u-H[tr]+l^ 

transforms the equation (jlip into the fixed-point problem 

u = -QoTz(^u- H[f:r] -non(^u- H[tr] + . (12) 

up to the CO- kernel error term in /C. The remaining task is to show that the mapping 

given by the right hand side of (I12p is a contraction mapping onto a neighbourhood of zero con- 
taining i7[Er] — ^. Finally, we will show that a well-behaved solution of equation (jlip up to an 
error term can be found for any choice of (r, W, H) where H is sufficiently small, W is contained in 
a compact subset of Hgev ^g-^ \ ^ where A is the subset of IlgGV ^<?-^ consisting of those W for 
which the unit vectors associated to any pair of neighbouring vertices point in opposite directions 
(after parallel transport along the geodesic segment connecting the vertices in question so that 
the vectors can be compared), and r is sufficiently small but all its components are positive and 
bounded away from zero by a fixed amount depending only on r. 

3.3 Structure of the Mean Curvature Operator 

In this section, we collect all the various structural facts about the mean curvature operator that 
will be used in subsequent sections of this paper. We can more or less quote exactly the analogous 
results from [1] since the setting here is in most respects identical to the setting of [Ij. Nevertheless, 
a brief development of the results will be presented here in order to give the reader a sense of the 
structure of the partial differential equation that must be solved in order to deform T,r into a CMC 
surface. 

Normal coordinate charts. Consider a geodesic normal coordinate chart centered at a point 
p € M. We can express the background metric in this chart as 

g:=g + P := [Sij + Pij{x)) dx* (g) dx^' 

where g is the Euclidean metric and P is the perturbation term. It is well known that 

M^) = ^ E^*'^-^^')^'^" + ^ E Ra.m-Mx'x'^x'' + Odixf ) . (13) 

l,m l,m,n 

Where R,,, := Rm(^, £, ^, and R,,,^ := V^Rm(^, £, ^, are components of 
the ambient Riemann curvature tensor Rm and its ambient covariant derivative VRm. 
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Surfaces in normal coordinate charts. Now let S be a surface (perhaps with boundary) 
contained in this chart. We can now use the above expansion to derive corresponding expansions in 
terms of the background curvature at p for any geometric object defined on S. Here and in the rest 
of the paper, let h, T, V, A, A^, B, H be the induced metric, Christoffel symbols, covariant derivative, 
Laplacian, unit normal vector, second fundamental form, and mean curvature of S with respect 
to the metric 5, and let h,T,'V,A,N,B,H be these same objects with respect to the Euclidean 
metric. Near a point x G T,, let {Ei,E2} be a local frame for TS induced by some coordinate 
system and denote by Y := gfj position vector. Define 



'Poj 
Vijo 



P{N,N) 
P{N,E,) 

\{{EiP){Ej,Et) + {EjP){Ei,Et) - (EtPm^Ej)) 
\{{E,P){Ej,N) + {E.jP){Ei,N) - {NP){Ei,Ej)) . 



Straightforward geometric calculations now give us the following results. 

Lemma 3. Let Ti be a surface belonging to a geodesic normal coordinate chart of M where the 
expansion ()13p is valid. Then the induced metric of S and the associated Christoffel symbols satisfy 



K,=K,+Vij and h'^ = -W'Wst + 0{\\Yf ) and T.^k = Tijk + V^Jk + VokB, 



the normal vector of S satisfies 

N = = (1 - hVoo)N - U^Vo,E, + Odlyf) , 

[l+Vo^-h^iVoiV^^r^ 

the second fundamental form of S satisfies 



B.J = (1 + Too - U^Vo.Vo,y'^B., + ^-^0 h^'W.^i 

{l + Voo-h^^iVoj) ^ 

= (1 + ^Voo) Bij + Vijo + Bij{Y, B, N, ^1, ^2) , 
and finally the mean curvature of S satisfies 

H={1 + \Voo)H + h'^Vijo - B'^Vij + n{Y, B, N, ^1, ^2) • 
In the last two sets of identities, Bij and Ti are functions with 

max\BijiY,B,N,Ei,E2)\ + \n(Y,B,N,Ei,E2)\ < C\\Yf{l + \\Y\\\\B\ 
for a constant C depending only on the curvature tensor of the ambient manifold. 
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Normal graphs in normal coordinate charts. Let us next suppose that the surface S in the 
geodesic normal coordinate chart where the expansion (jl3p is valid is deformed in the following 
way. Choose a function / : S ^ M and define /ij : S ^ R'^ to be the normal deformation of S 
by / with respect to the normal vector field of S calculated using the background metric g. Put 
Tif := We would like to obtain a workable expression for the mean curvature of T,f and 

understand its dependence on the function / and the expansion of the metric g. 

However, straightforward comparison is not possible because in order to express Sj- as a normal 
deformation of S with respect to the Euclidean metric, we can not use a different normal graphing 
function. Let /xj : S ^ to be the normal deformation of S by / with respect to the normal 
vector field of S calculated using the Euclidean metric. Put := //j-(S). The next lemma shows 
that we can replace by /i/o(S) for a new function /o : S ^ M and then relate |/|c2>« to \ fo\c<2.a. 

Lemma 4. Let Ti be a surface belonging to a geodesic normal coordinate chart of M where the 
expansion is valid and let := Under the assumption that diam{Ti) is sufficiently 

small and |/|||i5|| ^ 1 on S, then there exists a function /q : S — > M so that = fiff^{T,). Moreover, 
/o satisfies the estimate 

l/olc2."(Be) < C'l/lc2."(B2e) 

where C is a constant that depends only on the geometry of M in the coordinate chart, while Bg 
and b20 is any pair of concentric balls of radii 6 and 29 contained in the coordinate chart. 

We can now analyze the mean curvature operator and the second fundamental form of a surface 
of the form := in a geodesic normal coordinate chart, where f : T, ^ R. We start with 

the fact that the second fundamental form and mean curvature of S j with respect to the Euclidean 
metric can be decomposed as 

B[±f]=B + m{f) + B'-'\f) ^^^^ 
H[tf]=H + t{f) + Q{f) 

where t:=K+ \\Bf and ^(^) are linear operators while Q and B^ ' are second-order differential 
operators that are quadratic and higher in their arguments. Expressions for these operators are 
well known and are given in \l\ . Using Lemma [3] we can now derive from (jl4|) the analogous 
decomposition of the mean curvature of S/ with respect to the background metric g as well as 
obtain estimates for /. This is carried out in [IJ and the result is the following. We note that the 
requirement + ||V/|| <^ 1 appearing below is actually very natural and we will show later 

that it holds for the functions we will be considering. 

Lemma 5. Let Ti be a surface belonging to a geodesic normal coordinate chart of M where the 
expansion (jl3p is valid and let T,f := fifiTi). Then the mean curvature ofTij with respect to the 
background metric g can be expanded as 

H[±f]=H + C{f) + Qif)+n{f) 
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where £ is a linear operator, Q is a quadratic remainder and Ti is an error term. Furthermore, 
the following estimates are valid. First, C satisfies 

\C{u) - C{u)\ < C(l + ||y||||-B||)(|u| + ||y||||Vn|| + ||yf iiv^nii) . 

Under the assumption that + ||V/j|| <^ 1 fori = 1,2, then Q satisfies 

|Q(/i) - Q(/2)| < C\fi - f2\ -m^xil fi\\\Bf + \\Vfi\\\\Bf + ||V/,||||V^|| + || V^/^ll P||) 
+ C||V/i - V/2II • max {\fi\\\Bf + ||V/i|||iS|| + |/^|||VB|| + \\V^f^\\) 
+ Clival - V2/2II • max [m\\B\\ + \\V fi\\) 

+ C\\V^fi - V2/2II • max ||V/,f + C||V/i - V/2II • max || V/i || || V^/ill 

i i 

+ C\fi - /2I • max WVfiW + C7||V/i - V/2II • max \fi\ 

i i 

+ C{\fi - /2I + ||y||||V/i - V/2II) • max (l/il + ||y|||lV/i||) . 
Finally, under the assumption that + ||V/|| <^ 1 then Ti satisfies 

Mfi) - n{f2)\ < c\\Yfii + ||y||||^||)(|/i - /2I + ||y||||v/i - vmi + llvVi - V2/2II) . 

In all of the estimates above, C is a constant depending only on the curvature tensor of the ambient 
manifold at the center of the normal coordinate chart under consideration. 

3.4 Estimates of the Mean Curvature of the Approximate Solutions 

The first step in the proof that the equation //[Sr] = f can be solved up to an error term is to 
estimate -f^[Sr] — f o a ■ need only trivially modify the estimates found in [Ij of 

this quantity for use in the present setting. Therefore we give only a very abbreviated outline of 
the proof here. In the what follows, set e := maxje^ : q G V} U {e^b : q^ G V**} and recall that 
r, = 0(re3/4). 

Proposition 6. Suppose v G (1,2). The mean curvature ofT,,. satisfies the estimate 

2 



r 



< Cmaxjr^-^ r^-'^e^/^-^'^/'^j 



for some constant C independent of r and e. Moreover, this estimate is uniform in (r, VF,H) 
provided these satisfy the requirements set out in Section \3. 



Proof. The proof follows three steps. The first step is to derive pointwise estimates for the mean 
curvature and second fundamental form of S,. with respect to the Euclidean background metric 
in the spherical and neck regions in the geodesic normal coordinates in which they are defined. 
For this purpose we use the well-known explicit expressions of these quantities that are available. 
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The second step is to convert these estimates into pointwise estimates for the mean curvature with 
respect to the background metric g using Lemma H] and Lemma \E\ The third step is to deduce 
the estimate of the C^!_2 norm of H[tr] - ^. The result is as stated above, where the uniformity 
of the estimate in W follows from the first step since the Euclidean mean curvature and second 
fundamental form of the building blocks is uniform in W provided the unit vectors pointing to the 
neighbouring building blocks are not too close together. □ 

3.5 Analysis of the Linearized Mean Curvature Equation 

The second step in our proof is to find a parametrix TZ satisfying C oTZ = id -\- £ where £ maps 
into a subspace of functions /C. We construct this parametrix by patching together solutions of 
the Euclidean linearized mean curvature equation in each geodesic normal coordinate chart used 
to define the constituents of using a procedure that is in broad terms completely analogous to 
the construction found in p]. However, because we have not imposed any symmetries whatsoever 
on Sj., the approximate co- kernel is much larger now than in [T] and this makes the construction 
of our parametrix rather different in the details. Thus we present a more thorough proof. 

Proposition 7. Let v E (1,2). There is an operator TZ : C|^^2(^r) Cu'^'CEr) that satisfies 
CoTZ = id — £ where £ : C'^'^2{^r) ^ md iC is a subspace of functions on T,r that will be defined 
below. The estimates satisfied by TZ and £ are 

\TZ{f)y,. + \£{f)\cO,^ < C|/|^o,. 

for all f £ C'^'^2i'^r), where C is a constant independent of r and e. Moreover, this estimate is 
uniform in (r, VF, S) provided these satisfy the requirements set out in Section \3. 

Proof. Let / G C|^^2(^r-) be given. The task at hand is to solve the equation C{u) = f + £{f ) for a 
function u G C^'"(Er) and an error term £{f) G K,. To begin, introduce four radii cJi < (T2 < cjs < 
(74 ^ r with the property that the supports of the gradients of the cut-off functions x*^ and x*^ 
do not overlap for i / j. 

Step 1. Choose g'' G and let f^eckq' ■~ J'^'^eckq*'' "^^^^ function can be viewed as a function 
of compact support on the standard scaled catenoid re^bN . Now consider the equation Cmckiu) = 
fneck,q^ ou rs^iN , where Cmck is the linearized mean curvature operator of re^^N with respect to 
the Euclidean metric. Then the pull-back of C to re^bN is a small perturbation of Cneck- By the 
theory of the Laplace operator on asymptotically flat manifolds, the operator Cneck is surjective 
onto C^!"2 (^^g^ -^) when v G (1,2). Hence there is a solution u^^^^^^b as desired, satisfying the 
estimate \u^^^j^ gi'lc^-" ^ C'[/[j^o,c« where these norms can be taken as the pull-backs of the weighted 
Holder norms being used to measure functions on S^. Finally, extend u^^^f^ ^b to all of by the 

definition U^^^^^^b := U„^^k^^bX2ck,q^ ^^leck ■= E^^bgy' ^neck,q^- 
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Step 2. Choose g G V and suppose that the spherical building block Sq := is built from 
§2 \ {pi, ■ ■ ■ ,Png}- Let fsph,q ■= (/ - C,{uneck))x1lh,q- ^hcn fsph,q is a function of compact support 
on Sq and thus can be viewed as a function of compact support on the sphere \ {pi, . . . ,Png}- 
Now consider the equation Csph{u) = fsph.q-, where C^ph is the linearized mean curvature operator 
of S'^ with respect to the Euclidean metric. Then the pull-back of C to the sphere is a small 
perturbation of Csph- It is not a priori possible to solve the equation Csphiusph,q) = fsph,q on 
§^ because Cgph possesses the three-dimensional kernel /Cs2 := span{ Ji, J2, J3}. Let f^f^q be the 
projection of fsph,q to the orthogonal complement of /C§2 with respect to the Euclidean L^-inner 
product. Now there is a solution of the equation Csph{usph,q) = fjphq satisfying the estimate 
I^isp/i,g|c2."{§2) < C'|/sp/i,glco."(s2)- One thus has \usph,q\c'^,'^{rS.'^) < Ccjg"^ I /Ij^o^^ when scaled and 
pushed forward to Sq. 

A modified solution satisfying weighted estimates on Sq can be obtained by exploiting the 
behaviour of the Taylor series expansion of Usph,q near the points pi. We can write 

n<j 2 

Usph,q '■= Vsph,q + {^q,i + ^ ^ ^q,i^q,i^'^q,i 
i=l s=l 

where Vgph^q satisfies \vsp}h,q{x)\ ^ Ca2~'^d\si{x,piy'\w\^o,a near pi and the remaining quantities are 
as follows: real numbers aq^i := Usph,qiPi), real numbers 6^ j which are components of the vector 
yUsph,qiPi) in a fixed basis of Tp-S^, associated functions i^^ that are linear in the distance from 
Pi, and finally functions rjq^i that equal one near pi and vanish a small e-independent distance away 
frompj. Furthermore, Uq^i and b^^ satisfy ||a|| + ||6|| < C(T2~^|/|^o,a . Finally, by adding the correct 
linear combination of the Js to Ugph^q, we can always ensure that the a,b quantities associated to 
one of the points pi vanish. 

We now extend each Usph,q to S^, and combine them as follows. We let Usph '■= Vsph + ^ where 
^sph ■■= Eq&;X?ph,qVsph,q and A := Eq&;T,Zi + EsK/g,i)^i'i^neck,p,- The function u^ph is 
now defined on all of and we have the estimates \vsph\(j^<<^ + o'2~'^(ll^ll + ll^ll) — ^\f\c°-" • 

Step 3. Let u^^) := v,ph + Uneck and S^^\f) := - T.qevX7phj!ph,q " EgGV ^^ere 
fsphq projection of fsph,q to spanjjq^s : s = 1,2,3} with respect to the Euclidean L^-inner 

product and Jq^s is the push- forward of the function Jg to Sq. By collecting the estimates from 
Steps 1 and 2, one has lu*-^-* |(^2,q < C\f\^o,a . Then using the same arguments as in [Ij, we find that 

|/:(n«) - / - £^'\f)\cO,^^ < 0\f\cO,^^ and \£^'\f)\cO,^ < CarVlc"- ^^^^ 

where 9 can be made as small as desired by adjusting ai, . . . , (T4 and e suitably. The consequence 
is that one can iterate Steps 1 and 2 to construct sequences u*-"^ and £^^\f) that converge to 
u := TZ{f) and £{f) respectively, satisfying the desired bounds. □ 

The definition of the finite-dimensional image of the map £ : C^'^2(^r) —^iCisa by-product of 
Step 3 of the previous proof. 
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Definition 8. The approximate co-kernel of the operator C is the subspace 

K, := span{x^p\ ^ : g G V and s = 1, 2, 3} 

e span{x''gpf^^gCsph{vq,iXnlck,pJ '■ where Sg has necks at pi,... ,p„J 

e span{x'',pf^ gCsph{vq,iXllck,p/q,i) ■ s = 1,2 and g G V where S'^ has necks at pi, . . . ,pnj • 

3.6 The Non-Linear Estimate of the Mean Curvature Operator 

The next step in our proof is to find estimates for the norm of the quadratic remainder term 
Q and the error term 7i appearing in the expansion of the mean curvature operator. Once again, 
we need only trivially modify the analogous estimates found in [1] for use in the present setting. 
Therefore we give only a very abbreviated outline of the proofs here. Since the desired estimates 
come from combining Lemma U] and Lemma El we must first justify the assumption required there. 

Lemma 9. Pick x G T,r. Then x belongs to one of the normal coordinate charts used in the 
construction of T,j. where the second fundamental form with respect to the Euclidean metric is 
B{x). At this point, the estimate 

r{\\B{x)\\\f{x)\ + \\Vf{x)\\)<Cr^\fy^^. 

holds for all f G C^''^(Sr-); where C is a constant independent of r, e and 6. 

Hence it is true that r(|/|||i3|| + ||V/||) <^ 1 can be ensured by keeping |/|(^2,c« small enough. The 
following estimates are a consequence. 

Proposition 10. There exists M > so that if /i,/2 G C^'"(Sr) for v G (1,2) and satisfying 
|/i|^2,c« + \f2\c<2,c < M, then 

|Q(/i) - Q{f2)\cl^^^ < Cr''~'^\fi - f2\cl'" max{|/i|^2,a, |/2|^2,a} 
|W(/i) - W(/2)le0.c < CrVi - f2y,^ 

where C is a constant independent of r and e. Moreover, the estimates are uniform in (r, M^, H) 
provided these satisfy the requirements set out in Section \3. 

3.7 The Contraction Mapping Argument 

We are now in a position to solve the CMC equation up to a finite-dimensional error. Let E[r, e) := 
maxjr'^"'', 7'i-'^£3/2-3i//4| assume < e < ^ 1 (this will be justified a posteriori). The 
following estimates have been established. 

• The mean curvature satisfies — -l^o.a < CE[r,e). 

• There is a parametrix TZ satisfying ColZ= id — 8 where 8 maps into the finite-dimensional 
space it and \TZ{f)\^2,a + |f (/)|^o,a < C|/|^o,q for all / G C^!^2(^r)- 
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• The quadratic remainder satisfies |Q(/i) — Q(/2)|(70.« < Cr'^ — /2|(^2,q maxj { |/j|^2,a | 
for all /i,/2 S C^'"(Sr) with sufficiently small C^'"^ norm. 

• The error term satisfies \'H{fi) — TC{f2)\^o.a < Cr^|/i — f2\(f-,<^ for all /i, /2 G C^'"(Sr.) with 
sufficiently small C^'" norm. 

One can now assert the following. 

Proposition 11. There exists f := fr € (7^1*2(2,,.) and corresponding u := Ur ^ C^''^(l].r) defined 
by u := i7^(/ -H[tr\ + f ) so i/iai 

H[^jirf{tr)\-'^ = -8 (16) 

where £ (zK. The estimate |/|(^2,a < Cr^^E{r,e) holds, where the constant C is independent of r, 
e and is uniform in (r, W, H) provided these satisfy the requirements set out in Section \3.2l 

Proof. By the last three bullet points above, the map / i— > Mr{f) satisfies 

\K{fi) - K{f2)\clt\ ^ C{r^'^E{r,e)+r^)\fi - M^o^.^ 

where C is independent of r and e. Since r'^^^E{r, e) and r^ can be made as small as desired by a 
sufficiently small choice of r and e with < e < it is thus true that Mr is a contraction mapping 
on the ball of radius E{r, e) for such r and e. Hence a solution of p6|) satisfying the desired estimate 
can be found. The dependence of this solution on the parameters (r, W, H) is smooth as a natural 
consequence of the fixed-point process. □ 

4 Finding Exactly CMC Surfaces 

At this point, we have found a deformation of T,r{T,T,W,E) (with (r, VF, H) satisfying the usual 
conditions) into an almost-CMC surfaces, meaning that H[tr{T, r, W,E)] - ^ e )C where £ is the 
space defined in Definition [HI It remains to show that one can make a choice of F and (r, W, H) 
so that H[T,r(r,T,W,E)] — | = identically. The strategy for doing so is: first to relate the 
components of i?[Sr(r, r, W, H)] — ^ in ^ to the geometry of T and to the r and H parameters by 
means of the so-called balancing formula; then to use this formula to state conditions for which 
the equation H[tr{T,T,W,E)] - ^ = can be solved. A balanced surface is one for which these 
conditions hold. 

4.1 Derivation of the Balancing Conditions 

Suppose S := T,r(T,T,W,E) is a given initial surface and let Sj := ^^^(^^(r, r, VF, H)) be the 
almost-CMC surface generated by the function / := /^(t, VF,H) that solves the constant mean 
curvature equation up to a co-kernel error term belonging to /C. We now define two projection 
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operators. First, for each g G V, define tt^ : C^!^^2(^) ~^ ^ ^g(^) •= vrq^2(e), vrq^3(e)) 

where 



TTqAe) •= ^ ■ ^''sph,qJQ,s dVolg 



Second, for each g € V define functions tj ^ and ^ [, with the following properties: tj [, equals 1 

q g ,s g 

on the lb end of the neck at q" and is identically zero on the other end; while £^ ^ equal the linear 
functions and on the =b end of the neck at (in the coordinates used to define this neck) 
and is identically zero on the other end. Now define TT^b : C^;^2(^) ~^ 

V •= (v,o,+('^)' V,2,+('2)' V,o -(^)' V,i -(^)' V,2 -(e)) 

where 

V,s,±(e) - / e•^^P'^«W^^JdVol, . = 1,2. 

The following lemma gives the action of the projection operators TTq^g and vr^b ^ ^ on the basis 
for K. given in definition [HI It implies that if e € ^ and iTq{e) = TTqt,{e) = for all q and q^ &v'^ 
then e = 0. The proof is a straightforward computation. 

Lemma 12. The projection operators TTg^s o,nd vr^b ^ ^ satisfy the following properties. 

'^qAx7ph,g"J<l',s') = Sqg'Cr'^{6ss' + o(l)) 
'^g,s{Xsph,q'^sph{'nq',iXneck,pi)) = ^ 
{Xsph,q'^sph{'ng' ,iXneck,p/g' = ^ 

0{r~'^) is adjacent on the ± sirfe o/ g' 

otherwise 

-4^-2/ 



Cr (T-j^ ('^q^pi'^sO +o(l)) (7 is adjacent on the ± sirfe of q' 
otherwise 



'^q\s,±{Xsph,q'^sph{'nq' ,iXneck,pi)) 

U 

/ -X _ j C''^"'^0-f^((^gi,p.(5<js'+o(l)) q^ is adjacent on the ± side of q' 

''^ q\s,±\Xsph,q' '-'sphyVq' ,iXneck,p/g' ~ ^ 

I otherwise 

The next task is to use the expansions of the mean curvature and the estimates of the function 
/ to derive a formula relating 7r^,(i?[Sj] — 2/r) to the geometry of S. 

Proposition 13. The mean curvature ofTif satisfies the following formuliB . First, 

^ ^ adjoining ^ X TQj^ 
necks 

+ 0{er^) + 0{re^) + 0{r'^E{r,e)) 
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where the sum is taken over necks adjoining q and [Xgi,]s is the s- component of the unit vector at q 
pointing in the direction of the neck at q^ . Here Ci, C2 are numerical constants and R is the scalar 
curvature of M . Second, 

V,^,± (^P/l - 7) = E «f + ^"'^r'^(S,0 + 0{r-^e^) + OK-V-ii^(r, e)) (17b) 
^ ^ i=i 

where is a non-zero multiple of the i^^ component ofB^b, the function ^ is a function bounded 
uniformly by an r -independent constant times ||H^b|p, and the matrix Mis is given in (I19p below. 

Proof. The derivation of formula ()17ap is identical to the derivation given in [Ij . This derivation is 
based on the fact that the expansion of the mean curvature in terms of the background geometry 
given in equation [TS] reads 

H ={l + iRm(iV, y, iV, y) + 3^VyRm(iV,y,7V,y))F 

- (iRm(^i, y, Ej,Y) + iVYRm(S,, Y, E,,Y))B'^ 

„ _ _ o o (18) 

- |Ric(y,A^) - iVYRic(y,iv) + 3^v^Ric(y,y) - iv^Rm(iv,y,iv,y) 
+ o(||yf(i + ||y||p||)). 

which yields the leading terms in formula (|17ap plus error terms when integrated over T, f . 

We therefore proceed to the derivation of formula Il7b[ Choose q^ (z and consider the neck 
region A/J^^ Then '^qb^s,± (-^[^/] ~ f) integral over the transition region of M^^,^ and its 

neighbour on the it side, let us say that the neighbouring sphere is centered at g G V. Let p be 
the point in this sphere where the neck is attached. Using Lemma [Hand Lemma [5] as well as the 
definition of Sj, we can now write fl as a normal graph with respect to the Euclidean 

metric over an annulus A of radii 0{a\) about p. The graphing function is of the form fG \ A—^M. 
with G = G + ^-= + f where G is the generating function of the building block surface corresponding 
to q and $h is the deformation of the surface S n caused by the non-optimal matching to 

the neck region with parameters H^b , while / is the small deformation of / guaranteed by Lemma 
m The function is, to highest order in the magnitude of H, a linear combination of the form 
Yl^i=i'^t^i where $j := x^g^^ ^b5'(^«i and Xi is one of the translation, rotation and dilation 
vector fields, and af is a non-zero constant multiple of the parameter in H^b that corresponds 
to it. We define: Xi,X2,Xz generate translation along the directions, respectively, in 

the coordinates used to define the neck; X4 generates dilation; X^,Xq generate rotations in the 
x^, x^-directions, respectively. We also set ^ := t?^ for convenience. Hence 

+ r- 2^(H^b ) + 0{r-^e^) + 0{a\-^r-^ E{r, e)) 

where the function ^{■) is a quadratic remainder term, bounded uniformly by a constant indepen- 
dent of r and e. To conclude, we define the matrix Mis '■= f ^sph{^i)^sph{xZ^^. ^b^^b 5) dVolg and 
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times 1 or cos{9) or sm{9). Consequently 




where Cj, C'^ are constants and all other entries are zero. This is what we wanted to show. 



□ 



4.2 Balanced Almost-CMC Surfaces 

Suppose S := Y!,r{T,T,W,E) is once again a given initial surface and let := firfCEriT^TjW^E)) 
be the almost-CMC surface generated by the function / := /,.(t, W^S) that solves the constant 
mean curvature equation up to a co-kernel error term belonging to IC. The extra step needed for 
being able to deform S into an exactly CMC surface is to find a network of curves T and a value 
for the {t,W,E) parameters for which the projections 7rq^s{H[T,f] — 2/r) and vr^b — 2/r) 

can be shown to vanish for all g G V and g € V^ We will not be able to give a definitive answer to 
this question; rather we give a simple, testable set of conditions on T such that if these conditions 
hold, then the deformation to a CMC surface is possible. 

Balancing equations for the neck regions. We first consider the equations that must hold in 
the neck regions of S j and it turns out that these are easy to satisfy. There are two sets of equations 
for each v'^ corresponding to the it ends of the neck 7V^^ and these are given in equation 
()17bp . Now the quantity af is a constant multiple of the i*'' component of H^b. By the nature of 
translation, dilation and rotation of the catenoid, the functions have certain symmetries: 
^>5, <&6) corresponding to translation in the x^-direction and rotation in the x^- and x^-directions, 
are odd functions with respect to the center of the catenoid; whereas the functions <1>2, ^3, ^4, 
corresponding to translation in the x^- and x^-directions and dilation, are even functions with 
respect to the center of the catenoid. Consequently if we set af := ai for each i = 1, . . . , 6, then 

= —ai =02 '^3 = '^3 ^4 =04 ^i^d ^6 ~ • 

Therefore the two sets of equations pertaining to the neck A'^^^^ can be combined and now read 

6 




(20) 



i=l 
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where 



.-3^-2 



C[r-^<jf\og{l/ai) 



Mi, 



v 



C[r-^a^'\og{l/ai) 



Since Mis is an invertible matrix and ai = 0{r) we can solve equation (j20p for (ai,...,a6) 
near (0, . . . , 0) when r is sufficiently small. The solution depends parametrically on the (r, W)- 
parameters of the initial surface Ilr(r, t, W, H). We will continue to denote by Sj the deformed 
initial surface with the choice of H determined here. Moreover, one can check that the parameter 
values H^b that follow from our solution of ([20]) generate an o(r) deformation of Tj^ir^ W,Q) in the 
C^'" norm. Let us denote these parameter values collectively by H* := H*(r, r, W) and the surface 
we get as Sr(r, r, W, H*). 



Balancing equations for the spherical constituents. We now turn to the equations that 
must hold in each spherical constituent of Sj. Quoting equation (jl7p . we know that we must find 
(r, W) so that 

a o 

0= ^ Cire^,[7Vg.],-C2r4 — (g) + (!?(er3) + 0(re2) + 0(r2i?(r,e)) VgEV (21) 

adjoining 
necks 

where the sum is taken over necks adjoining q and [A^^bjs is the s-component of the unit vector 
at q pointing in the direction of the neck at g''. At this point, we make use of the existence 
conditions that we have imposed on F and explained in the introduction. We have shown that 
the leading-order terms in (j2ip vanish when the curves in T and the choice of neck sizes satisfy 
Conditions (2) - (3). By Condition (1) that the error term involved in approximating the leading- 
order part by V-y7 — Qr^VR o 7 is as small or smaller than the 0{er^) + O(re^) -|- 0{r'^E{r,e)) 
terms already present above. Here VL := C2/C1. We obtain a discrete family of balanced surfaces 
(r*, r*, 0, S*(r*, r*, 0)) because only for r belonging to a discrete family of shrinking intervals 
converging to zero is it possible to place an integer number of spheres and necks along each curve 
in r to within an error smaller than the leading order terms of (j2ip . Finally, the non-degeneracy 
condition, Condition (4), guarantees that the map taking the remaining VK-parameters onto the 
image space of all the vr^^s projections is locally surjective. Hence for sufficiently small r, we can 
find W := W* near zero so that the right hand side of ()2ip vanishes identically. Our desired CMC 
surface is S^(r*, r*, VF*, H*(r*, r*, VF*)). □ 
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